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; [ We consider an invariant formulation of the system of Maxwell's equa- 

tG • tions for an anisotropic medium on a compact orientable Riemannian 3- 

manifold (M, g) with nonempty boundary. The system can be completed 
to a Dirac type first order system on the manifold. We show that the 
Betti numbers of the manifold can be recovered from the dynamical re- 
sponse operator for the Dirac system given on a part of the boundary. 
In the case of the original physical Maxwell system, assuming that the 
entire boundary is known, all Betti numbers of the manifold can also be 
determined from the dynamical response operator given on a part of the 
I/"") ■ boundary. Physically, this operator maps the tangential component of the 

electric field into the tangential component of the magnetic field on the 
l/~) ■ boundary. 
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1 Introduction 



Recently there has been a lot of interest in inverse problems for Maxwell's 
equations in Euclidean domains in IR 3 and on compact Riemannian manifolds, 
see [i 13 HI EH [2U [23] . In a smooth bounded domain M C M 3 , Maxwell's 
equations are given by 

cuT\E(x,t)=-B t (x,t), 

curl H(x,t) =D t (x,t), 
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where E and H are the electric and magnetic fields, and B and D are the 
magnetic flux density and the electric displacement. The fields E and D, and 
similarly, the fields H and B are related by the constitutive relations, 

D(x, t) = e(x)E(x, t), B(x, t) = n(x)H(x, t), (1.2) 

where the electric permittivity e(x) and the magnetic permeability fJ.(x) are 
C°°-smooth positive-definite 3 x 3- matrix valued functions on M. The initial 
boundary value problem for the time dependent Maxwell's equations consists 
of (|l.ip . (|1.2p together with the conditions 

E(x,t)\ t =-r f = 0, H(x,t)\ t =-r f = 0, 

(1.3) 
n x E\ dM xR- = f, 

where n is the unit exterior normal to dM, and Tf > is such that f(x, t) = 
for t < —Tj. The inverse problem associated with (jl.ip . (|1.2p . and (jl.3p . is the 
problem of reconstruction of electromagnetic parameters e(:r) and /u(x) from 
the knowledge of the response operator 

R:nxE\ dMxR _^nxH\ dMxK _. (1.4) 

From the point of view of modern electrodynamics and classical field theories, 
it is natural to adopt an invariant approach to Maxwell's equations, where the 
domain M is replaced by a general 3-dimensional smooth compact oriented con- 
nected Riemannian manifold, and the vector fields E, H, D, and B are viewed as 
differential forms, see [26J. The geometric inverse problem is then to determine 
the unknown manifold M, together with the electromagnetic parameters, from 
the response operator (|1.4p . which is now defined in terms of boundary traces 
of the corresponding differential forms. See also |17[ ITS] , where the problem of 
the reconstruction of a Riemannian manifold from the Dirichlet-to-Neumann 
operator for harmonic functions, has been studied. 

In the context of time-harmonic Maxwell's equations in an isotropic setting, i.e. 
when the parameters e(x) and fJ,(x) are scalar, the inverse problem for bounded 
domains in M 3 was solved in |21j , see also (BJ [19} [23] . Much less is known in 
the anisotropic case. To the best of our knowledge, the positive results in this 
direction have only been established in the case of an anisotropic medium of a 
special type, characterized by the polarization independent velocity of the wave 
propagation. In terms of the electromagnetic parameters, this amounts to the 
existence of a(x) > such that e(x) = a(x)fi(x). In this case, under a certain 
geometric condition, it is shown in [14] that, if the conformal class of e(x) 
and n(x) is known, the stationary boundary measurements identify uniquely 
the conformal factors. There are also counterexamples for uniqueness of time- 
harmonic inverse problems involving very anisotropic and degenerate material 
parameters [11[ 112] , In [16], the inverse problem for Maxwell's equations in 
the time domain for an anisotropic medium was studied, still assuming that 
the wave propagation is independent of the polarization. It was shown that 
the Riemannian manifold and the electromagnetic parameters can be recovered 



from the dynamical response operator similar to (jl.4p . given on a finite time 
interval. See also [I] for reconstruction of the wave speed. 

In this paper we shall be concerned with the case of a general anisotropic 
medium. Specifically, working in the geometric setting of Maxwell's equations 
on a manifold M, we are able to recover the Betti numbers of the manifold 
from the dynamical response operator, given on an open subset of the bound- 
ary. This can be viewed as the first step in attempting to reconstruct the 
geometry and topology of the underlying manifold, in the full generality of 
the anisotropic case. Let us remark that in the isotropic case, as well as in 
the case when e{x) = a(x)fj,(x), a(x) > 0, the reconstruction of the manifold 
and the electromagnetic parameters is based on controllability results, which in 
turn rely crucially on generalizations of the Tataru unique continuation theorem 
[9l [16]. In our opinion, the main obstacle in the study of the inverse problem 
for the general anisotropic Maxwell system is due to the fact that such unique 
continuation results do not seem to be available in this case. 

We would like also to mention the paper [2], where the reconstruction of the 
Betti numbers of a manifold from the Dirichlet-to-Neumann operator for the 
Hodge Laplacian on differential forms is studied. 

The plan of the paper is as follows. Section 2 is devoted to the description of 
our geometric setup, including the completion of the Maxwell system to a Dirac 
type elliptic system, and contains the statement of the main results. We also 
discuss examples that illustrate the significance of our results for the determi- 
nation of the topological structure of an unknown object from the boundary 
measurements. In Section 3, we prove the identifiability of the Betti numbers 
in the complete Maxwell case, while in Section 4 we establish our results for 
the physical Maxwell system. 



2 Preliminaries and statement of the main results 

2.1 Invariant definition of Maxwell's equations 

Let (M, <7o) be a smooth compact oriented connected Riemannian 3-manifold 
M with dM ^ 0. We shall first rewrite equations (jl.ip . (|1.2|) . in the anisotropic 
case, using the language of differential forms . In doing so, we shall follow 
closely [16], where the case e(x) = a(x)fi(x), a(x) > 0, is considered. 

Let A T*M, k = 0, 1, . . . , 3, be the bundle of the k-th exterior differential forms 
and AT*M be the full bundle of differential forms. Denote by C°°(M, A k T*M) 
the space of smooth real exterior differential forms of degree k. 

Define the fiberwise duality between 1-forms and vector fields, 

b : C°°[M,TM) -> C°°(M, A 1 T*M), X\Y) = g (X,Y), 
or in a coordinate system for X = a l -^-j, X = go^a^dx 1 . This map is bijective 



and has the following properties |25j : 

(curlX) b = * dX\ (dwXf = * d* X b , 

where 

d : C°°(M,A k T*M) -> C°°(M, A fc+1 T*M) 

is the exterior differential and *o is the Hodge operator with respect to the 
metric go, acting fiberwise, 

*o : C°°(M,A k T*M) -> C°°{M, A 3 ~ k T*M). 

We define the 1-forms £ = E° and % = H and the 2-forms B = *o-B and 
T> = *qD . Using the identity *o*o = id, valid in the 3-dimensional case, we 
can write Maxwell's equations (jl.ip in terms of differential forms as 

d£ = -d t B, dU = d t V. (2.1) 



Consider now the constitutive relations (jl.2p . We shall determine a metric g e 
such that the Hodge operator with respect to this metric, denoted by * e , satisfies 



V = * (eEy = * e £. (2.2) 

a 



In local coordinates {x l ,x 2 ,x 3 ), we have eE = e\E k -£j, (eEf = go t ij€ k E k dx % 



and thus, the middle term of (|2.2p yields 

* {eE) b = * (g ,i j e{E k dx i ) = -y/det{g )g^g Q , ij € j k E k si pq dx v A dx q 



= -^det(g )e{E k s jpq dx p A dx q , 

where s\ pq is the Levi-Civita permutation symbol. The right hand side of (J2.2 
implies that 

* e £ = * e (g 0)ik E k dx i ) = -Vdet(g e )gi J g ,ikE k s jpq dx p A dx q . 

Hence, (|2.2p is valid if we set 



\/det(ge)gl J go,ik = y /det (5'o)4- 



By taking the determinants of both sides, we get 



\/det(p £ ) = det(e) Vdet^o). 
Defining 

Jj — J n ki 

9 * ~ det(e) £k90 ' 
we see that (|2.2p is valid. Similarly, we see that for the metric tensor grf 

aife^S*' we have 

B = *o(vHf = ^H. 



Hence, the constitutive relations take the form 

V(x,t) = * e S(x,t), B(x,t) = * ll %{x,t). 
We consider the waves that satisfy the initial conitions 

B(x,t)\t=- T = 0, D(x,t)\ t= _ T = 0, r>0, 
Applying the divergence operator to (|1.1|) . we have 

div B(x,t) = 0, div D(x,t) = 0, tel, x € M. 
In terms of differential forms these equations imply that 

dB = 0, dV = 0. (2.3) 

In the further considerations, we will use only the pair (£, B) and denote it by 
(w 1 ,^ 2 ), where w 1 = £ and w 2 = B. The compatibility conditions (|2.3p imply 
that 

du 2 = 0, d* e uj 1 = 0. (2.4) 

It follows from (J23J that 

Let us consider the following codifferentials, 

8 e ^uj 2 = * e d * l _ l u 2 , 5^ e uj k = - * ll d* e uj k ,k = 1,3. (2.6) 

Then Q27JJ) and (g3]) yield 

u] = 5 € ^uj 2 , 5^ e uj l = 0, 
co t = —dui , du = 0. 



2 — — dcj , LJ t [ = * e d* fJ- uj 2 . (2.5) 



(2.7) 



These equations are called Maxwell's equations for forms in the divergence free 
case on a Riemannian manifold M. 

We shall now extend the above equations to the full bundle of exterior differ- 
ential forms AT*M. To this end, we introduce auxiliary forms, w° £ C°°{M) 
and w 3 G C°°(M, A 3 T*M), which vanish in the electromagnetic theory, by 

UJ t = S^efjJ , U) t = —du . 

Since w° = and w 3 = in the electromagnetic theory, we can modify equations 
(12771) to have 

uj t = —dco + 5 e ^uj , w t = —dui , 
w 2 = -doj 1 + <5 M , e w 3 , w t ° = S^eU 1 , 

or, in the matrix form, 

uj t + Duj = 0, (2.8) 



where uj = (u° , u 1 , u 2 , w 3 ) and the operator D is given by 

/0 -<W \ 

d -Ve ' ^ j 

\0 d / 

Equations (12.80 . (12.9(> are called i/ie complete Maxwell system. Notice that the 
operator D is of the Dirac type. 

2.2 Function spaces 

Define the L 2 -inner product in the space C°°(M, A k T*M) as follows, 

.k „fc^ / , .k a , „fe 

r2 ~~ 



K.TTjia = / w*Avj*, A: = 0,2, 



(u>W) L? = I u*A* e Ti*, fe = 1,3, 

•/AT 

and denote by L 2 (M,A k T*M) the completion of C°°(M,A k T*M) in the corre- 
sponding norm. In the complexified case, we take the corresponding sesquilin- 
ear extension of the inner product. We denote by H s (M,A k T*M) the standard 
Sobolev space of k- forms. 

The natural domain of the exterior differential d in L 2 {M,A k T*M) is 

H(d,A k T*M) = {uo k G L 2 (M,A k T*M) : duj k G L 2 (M, A k+1 T*M)}, 
and we define 

H(S e ^,A k T*M) = {uo k G L 2 (M,A k T*M) : <5 £iAt w fc G L 2 (M, A fc ~ 1 T*M)}, 
and similarly for <L e . 

Let z* : C°°(M, A k T*M) -»■ C co (dM,A k T*M) be the pull-back of the imbed- 
ding i : <9M —)■ M. Then we define the tangential trace of /c-forms as 

t :C°°(M,A fc T*M) ^ C°°(dM,A k T*M), tu k = i*oo k , k = 0,1,2, 

and the normal trace as 

n : C°°(M,A k T*M) -+ C°°(9M, A 3 ~ k T*M), 



nU } k = i*U e u k ), k = 1,3, nw 2 =i*(* M w 2 ) 



Set 



{tcu k ,m] k+1 }= [ tco k Anr) k+1 , k = 0,1,2. 
With this notation, Stokes' formulae for differential forms can be written as 

(doj^r, 2 )^ - (w 1 ,^ 2 )^ = (tw 1 ,^ 2 ), (2.10) 

[dw 2 ,rf) Ll - (w 2 ,<5 Mie r/ 3 ) L 2 = (tw 2 ,n77 3 ). 



Using d2HD and EH we get 

(Du, f]) L 2 + (w, Drj) L 2 = (tu, nry) + (try, nw), (2.11) 

where tu = (tea , tea 1 , tea 2 ) , nw = (nw 1 , nw 2 , nw 3 ), and 

(tw,nry) = (to; ,™? 1 ) + (tw\nry 2 ) + (tw 2 ,nry 3 ). 

Here we take w, ry E Ti, where 

H = H(d,A°T*M) x [H(d,A 1 T*M) r\H(5^A 1 T*M)} 
x [H(d,A 2 T*M)nH(S e ^,A 2 T*M)\ x H(5^ A 3 T*M). 

It will be convenient to write 5 to stand for both 5^^ and 5 e ^, when no risk of 
ambiguity is possible. There are well defined extensions of the boundary trace 
operators t and n to the spaces H(d,A k T*M) and H(5,A k T*M), see 



Lemma 2.1. The operators t and n can be extended to continuous surjective 
maps 

t : H(d, A k T*M) ->• H~ 1/2 (d, dM, A k T*M), (2.12) 

n : H(S, A k+1 T*M) -> H- l l 2 (d, dM, A 2 ~ k T*M), (2.13) 

where H- l / 2 (d,dM,A k T*M) is given by 

{u k e H- 1 ' 2 {dM,A k T*M) : du k G H~ 1 / 2 {dM, A k+l T*M)}. 

Let H t (d,A k T*M) stand for the kernel of (I2TT2D . and H n (5,A k+1 T*M) will 
denote the kernel of the operator (|2.13p . 

Using (|2.10p . we can verify the following result in a standard way, see also |16l 
Lemma 1.3]. 

Lemma 2.2. The Hilbert space adjoint of 

d : L 2 (M,A°T*M) -> L 2 (M, A X T*M), 

equipped with the domain Ht(d,A°T*M), is the operator <L g with the domain 
H(5 l j ije ,A 1 T*M). The Hilbert space adjoint of 

S^e : L 2 (M,A l T*M) -»■ L 2 (M,A°T*M), 

equipped with the domain H(5 il ^,A l T*M), is the operator d with the domain 
H t (d,A°T*M). 

It is clear that analogous statements hold for the operators d and 5, acting on 
forms of higher degree. 

We shall need the following result. 



Proposition 2.3. (i) The operator D, given by (I2.9P and equipped with the 
domain 

V{D) =H t (d,A°T*M) x [H t {d,A l T*M)nH{5^ e ,A 1 T*M)} 

x [H t (d,A 2 T*M)nH(5 € ^,A 2 T*M)} x H(S^ e ,A 3 T*M), 

is skew-adjoint on L 2 . 
(ii) The spectrum of the operator D with the domain T>{D) is discrete. 
(Hi) The operator D is an elliptic differential operator in the interior of M . 
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Proof, (i). Using the definition of the domain of the adjoint and Lemma 12.21 
we obtain that T>(D*) = T>{D). The skew-adjointness of D then follows from 
(|2"XL|) . which holds for uj,v £ V(D). 

(ii). In view of Gaffney's inequality \Zf>\ Corollary 2.1.6], 

H t (d,A k T*M) n H(5,A k T*M) = {u k e H l (M,A k T*M) : toj k = 0},k = 1,2, 

together with the Sobolev embedding, we conclude that the imbedding T>(D) ^-> 
L 2 is compact. Hence, the spectrum of T> is discrete. 

(hi). It suffices to show the ellipticity of D 2 . Since 5^ <e 5 t ^ = and <5 ei/i <5 M)(E = 0, 
we get 



D^ 



The operator D 2 enjoys the following coercive estimate, 

(D uj,u}) L 2(q M ) > Ci||u>|lHi(fiM) -C2|M| L 2 (nM ),Ci > 0, (2.14) 

where u = (wV 1 ," 2 ,^) and u k £ C^(M,A k T*M), k = 0,1,2,3. 
When proving (I2TT41 . notice that, for u 1 G C§°(M, A l T*M), 

((dd~^ t + Se^d)^ 1 ,u) 1 ) I ^ = H^eu; 1 !!^ + lldw 1 !!^. 

An application of Gaffney's inequality gives that 

ll^l^i < C(M)(|| W 1 || L2 + \\du l \\ L 2 + WS^co 1 ^) 

where C(M) > is a constant. The estimate (|2.14p follows, since the treatment 
of forms of degrees different from 1 is analogous. See also [7] for a different proof 
of coercivity. The ellipticity of D 2 now follows from the coercivity estimate 
(|2T4jl . see e.g. [20]. 

D 



2.3 Betti numbers and the Euler characteristic of a manifold 
with boundary 

Let (M, <7o) be an orientable compact Riemannian manifold of dimension 3 with 
boundary. The space 

% k (M) = {uj G L 2 (M, A k T*M) : du = 0, d * go oj = 0} 

is called the space of harmonic fields. Notice that this space is infinite dimen- 
sional for 1 < k < 2, see [25, Theorem 3.4.2]. Moreover, it is well-known that 
harmonic fields are C°°-smooth in the interior of M. The following two finite 
dimensional subspaces are distinguished in T-L k (M): 

U k D (M) = {uo G U k (M) : tw = 0} and 
U k N {M) = {to G H k (M) : i*(* go u) = 0}, 

which are called the Dirichlet and Neumann harmonic fileds, respectively. It 
follows from the Hodge theory that the dimensions of the spaces Tl k - ) {M) and 
% k N {M) are independent of the choice of the metric go- For our purposes, we 
shall have to specify the choice of the Hodge star operator in the definition of 
Ti k (M), according to the definition of the co differential given in (|2.6p . 

U 2 {M) = {to G L 2 (M, A 2 T*M) : du = 0, d * M u = 0}, 
H k (M) = {w G L 2 (M, A k T*M) : duj = 0, d * e u = 0}, k = 1, 3. 

Recall [10] that the space T~L%(M) is isomorphic to the &th homology group of 
the manifold H^M; R) and Tl k ^{M) is isomorphic to the fcth relative homology 
group Hk(M, dM;M.). The Poincare-Lefschetz duality states the existence of 
the following isomorphism, 

H k (M; R) ~ H 3 _ k (M, dM; M), k = 0, 1, 2, 3. 

The fcth absolute Betti number of the manifold M is given by 

p k {M) = dim H%(M), k = 0,1,2,3, 

and the fcth relative Betti number of M is defined by 

Pk(M,dM)=dimHn(M), k = 0,1,2,3. 

Being one of the simplest topological invariants, the Betti numbers carry a 
basic amount of information about the topology of a manifold in question. The 
Betti numbers /3q(M) and f3^(M) admit a particularly straightforward geometric 
interpretation. Namely, (3q{M) counts the number of the connected components 
of M and /^(M) gives the number of the oriented components of M without 
boundary. Assuming that the manifold M is connected, we have /3q(M) = 1 
and l3s(M) = 0. As for the first Betti number J3\{M), it is at least as large as 
the total number of handles of dM, see [5J Theorem 5.1.9]. 
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The Euler characteristic is defined by 

X (M) = /3 3 (M) - /3 2 (M) + ^{M) - /3 (M). 

It is known (SJ Corollary 8.8] that the Euler characteristics of a compact 3- 
manifold and its boundary are related by 

xipM) = 2 X (M). (2.15) 

Notice finally that if M is a connected compact orientable 3-manifold with van- 
ishing Euler characteristic, then either the manifold M is closed or its boundary 
is a disjoint union of tori. 

2.4 Boundary data for inverse problems 

Let r C dM be an open subset of the boundary dM. Consider the following 
initial boundary value problem, 

(d t + D)uj(x, t) = in M x R, 

tw| S MxM = / G CS°(R_,Cnr,AT*M)), (2-16) 

w|t=- T/ = 0, 

where 77 > is such that inf supp (/) > — 77. Following [16], we shall define 
a solution of (J2.16P in the following way. Let E be a right inverse to the trace 
mapping t such that Ef(—Tf) = 0. We set 

J{t) = Ef(t) - I e^-^ids + D)Ef{s)ds. (2.17) 

J-Tf 

Here e~ is the unitary group, generated by the self-adjoint operator D/i. 
Associated to the problem (|2.16|) is the response operator, 

R r : f i-> nu/|rxiL_- 

The first main result of this work is the following theorem. 

Theorem 2.4. Assume that we are given an open subset T C dM and the re- 
sponse operator Rr for any f £ C^°(1R_, C^°(T, AT*M)). These data determine 
the Betti numbers of the manifold M . 



Let us now return to the physical Maxwell's equations 

u] = 5 ttfl UJ 2 , d^eUJ 1 = 0, 
uj+ = —du] , duj = 0, 



(2-18) 

to; 1 = h e c^iR-^^ir^^M)), 
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As explained in [16], the solution to (|2.18p is obtained from (|2.17p by choosing 
the boundary source / in f|2. 16|) as 

f = (0,h,- I dh(t')dt'). 
The response operator for (J2.18P is defined by 

where u h is the solution to ()2.18|) . Notice that in the classical terminology 
of electric and magnetic fields, the response operator i?r maps the tangential 
component of the electric field n x .E|rxR_ to the tangential component of the 
magnetic field n x iJ|rxK_- 

Theorem 2.5. Given an open subset T C dM and the response operator Rr 
for any h £ C£°(M_, C^°{T, A 1 T*M)), the first absolute Betti number /3i(M) of 
the manifold M can be determined. 

Corollary 2.6. The knowledge of the boundary dM and the response operator 
R r , T C dM, for any h € C£°(R_,C§°(r, k x T*M)), determines the first and 
the second absolute Betti numbers /?i(M) and ^(M) of M . 

Corollary 12.61 follows from Theorem 12.51 together with (|2.15p . 



2.5 Examples 

The following two examples illustrate the significance of our results for the deter- 
mination of the topological structure of an unknown object from the boundary 
measurements. This may have applications to practical situations, where the 
structure of complicated voids in an unknown object is to be recovered. 

Example 2.1. Let M C K 3 be obtained from a large ball by removing a finite 
number of pairwise disjoint solid tori. Then the first absolute Betti number 
of M is equal to the number of the removed solid tori. Thus, measuring the 
response operator on a portion of the boundary sphere, we can recover the total 
number of the removed tori. 

Example 2.2. Consider a solid torus ST = S 1 x D 2 C M. 3 , where S 1 is a unit 
circle and D 2 is a closed two-dimensional disc. The boundary of ST is a two 
dimensional torus and since D 2 is contractible, it follows that the first absolute 
Betti number of ST is equal to 1. Let M be the connected sum of k copies 
of solid tori ST. Here we may recall that a connected sum of two manifolds, 
possibly with boundary, is a manifold formed by deleting a ball in the interior 
of each of the manifolds and gluing together the resulting boundary spheres. 
The boundary of M is a disjoint union of k copies of two-dimensional tori. It is 
known that for manifolds of dimension three and higher, the first absolute Betti 
number of the connected sum is the sum of the first absolute Betti numbers of 
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the summands. Therefore, the first absolute Betti number of M is equal to k. 
It follows from our results that performing measurements on a portion of the 
boundary of the manifold M, we are able to recover the total number of the 
solid tori. 



3 Proof of Theorem El 

3.1 Inner products 

Let w-'(i) = ojf(x,t) be the solution to (|2.16p . We shall need the following 
Blagovestchenskii type result, see [3j for such results for one-dimensional inverse 
problems. 

Theorem 3.1. For any f,h€ C^(M.-,C^{T,AT*M)), the knowledge ofTc 
dM and the response operator R? allows us to evaluate the inner products 

(u) f ' k (t),u h ' k (s)) L 2, k = 0,1,2,3, for s,t>0. (3.1) 

Proof. From (|2.17p . we obtain that 

cj f (t) = u> ft (-l), t>0, 

where f t = f(- + t + l), f t G Cg°(R_,Cg°(r,AT*M)). Therefore, the knowledge 
of the operator Rr is equivalent to the knowledge of the operator / — > nu/|rxR- 
To prove this theorem we also need the following fact, 

t{du k ) = d(tu k ), k = 0,1,2,3, 

see |25l Proposition 1.2.6]. This implies that 

n5 ejfJi LO =t(* e * e d * M u ) = dt(*uL0 ) = dmuj , 
nd^oj 3 =t(* M (— 1) * fl d* e w 3 ) = -dt(* e uj 3 ) = -dncu 3 . 

Set I k (s, t) = {u]f< k (t),u h > k {s)) L 2 , k = 0, . . . , 3. Then using Stokes' formulae, 
we get 

(d 2 -d?)i°(s,t) = (ojf'°(t),d^u h '°(s)) - (a t V'°(i),uA°(s)) 

= -(u;f> (t),5^duj h ' (s)) + (5^diuf> (t),uj h ' (s)) 
= (tu f '°(t),n(kj h ' (s)) - (tu h '°(s),nduj f ' (t)) 
= -(tu/^AnuA 1 ^)) + (tu; f ' (t),nS e ^w h ' 2 ( S )) 
+ (tco h ' (s),d t nio^(t)) - (tco h ' (s),n6 e ^ 2 (t)) 
= -(tw / '°(t),a s na; M (s)) + (tuj f '°(t),dnuj h ' 2 (s)) 
+ (tuA^dtnu/- 1 ^)) - (tw M (s), W- 2 (i)). 
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Similarly, 

(d 2 - d 2 )l\s,t) = {ujf\t),d 2 s uj h '\s)) - {%J\t),u h >\a)) 

= -{d s tuj h '°(s),nu; f > l (t)) - (tuj f ' 1 {t),d s nw h ' 2 (s) + dm; h ' 3 (s)) 
+ (d t tuj f ' (t),m; h ' 1 {s)) + (toj h ' 1 {s),d t noj f ' 2 {t) + dnoj f ' 3 {t)), 

(d 2 - d 2 )I 2 (s,t) = (^ 2 (t),d 2 ^ 2 (s)) - {dW>\t),^ 2 {s)) 

= -(d^ 1 ^) + dtu h '°(s),nujf> 2 (t)) - (tu f > 2 (t),d s nuj h > 3 (s)) 
+ (dttcj f '\t) +dtoj f '°(t),nuj h ' 2 (s)} + (toj h ' 2 (s),d t nLo f ' 3 (t)), 

(d 2 - d 2 )I 3 (s,t) = (^ 3 (t),d 2 co h ' 3 (s)) - (d^it),^ 3 ^)) 
= -(d s toj h ' 2 (s)+dtoj h ' 1 (s),noj f ' 3 (t)) 
+ (d t toj f ' 2 (t) +dtu J f ' 1 (t),noj h ' 3 (s)). 

Hence 1 (s,t), k = 0,1,2,3, satisfies an inhomogeneous one-dimensional wave 
equation in the unbounded region {(s,t) € M 2 : s > —Th,t > — Tf}, whose right 
hand side is determined from the knowledge of T and Rr . Since 

I k (-T h , t) = I k ( S , -Tf) = 0, d s I k (-T h , t) = d t I k ( S , -Tf) = 0, 

we can determine I k (s,t) in the entire region s > —Tf, t > — Tf. The result 
follows. □ 

3.2 Controllability result 

In the isotropic setting and the case when e(x) = a(x)/j,(x), a(x) > 0, one 
can use a generalization of Tataru's unique continuation theorem [9j [TB] to 
obtain controllability results with sources supported in a finite time interval. 
As already mentioned in the introduction, such unique continuation results do 
not seem to be available in the general anisotropic setting. Nevertheless, we 
shall next show that partial controllability results in the general anisotropic 
setting on an infinite time interval can be obtained using a unique continuation 
principle for elliptic systems. As shown below, this turns out to be sufficient 
for the reconstruction of the Betti numbers. 

Let T-Ld(M) := (B 3 , =0 ri k j(M) be the space of all Dirichlet harmonic fields, and 
let IT : L 2 (M,AT*M) ->■ Ud{M) be the orthogonal projection. 

Theorem 3.2. We have 

{nV(0)) : / G Cg (M_,CS (r,AT*M))} = ri D {M). 

Proof. Let rj € Hd (M) . If we prove that the orthogonality condition 
(w / (0), -r}) L i = for all / G C °°(IR_, C °°(r, AT*M)) 
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implies that 77 = 0, then the space {lW(0)) : / G Cg°(M_, C$°(T, AT*M))} is 
dense in Hrj(M). Since the latter space is finite dimensional, the claim follows. 

As Drj = 0, we shall view r/(x) as the solution to the following problem, dual 
to (I236|) . 

(-d t - D)u = 0, in M x R, 

tu\dMxR = 0, (3.2) 

u\t=o = r]. 

Using (J2.11J) . we have 

d t {uj f ,u) L 2 = -(Du f ,u) L 2 - (u f ,Du) L 2 = -(f,nu). 

Thus, 

o 

(f,nu)dt = -(u f (0), Tj)& + (u f (-T f ),u(-T f )) L2 = 0. 

- T f 

The choice of —Tf implies that 

(f,nu}dt = 



for all / g CS°(R-,C$°(T,AT*M)). Thus, nu = on T x R_. 

Now if r coincides with the whole boundary of the manifold M, then we are 
done, since H k D (M) n H%(M) = {0}, see (25J p. 130]. 

In the case when T is a proper open subset of dM, we proceed as follows. 
Notice that rj(x) solves the second order elliptic system D 2 rj = on M with 
zero Cauchy data on T, 

(trj, nr], t5r], ndrj), 

where 5r] k = 5^ >e r] k , k = 1,3, and 5r] k = 5 €:IJi r] k , k = 2. Thus, by the unique 
continuation principle for second order elliptic systems with diagonal principal 
part, see [13, Theorem 4.3], we get r\ = in M. 

D 

Corollary 3.3. Let U k : L 2 (M,A k T*M) -> U k D {M) be the orthogonal projec- 
tion onto the space of the Dirichlet harmonic k- fields. Then 

{n fe (u/' fc (o)) : / g c oo (R_,c oo (r,Ar , M))} = n k D (M), k = 0,1,2,3. 



3.3 Determination of the Betti numbers of the manifold 

Lemma 3.4. Let f,h G C£° (R_,C£°(T, AT*M)). Then given the response 
operator R^, it is possible to find the inner products 



(U k u f ' k (0),u h ' k (0)) L 2, k = 0,1,2,3. 
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Proof. Using (|2.10p , we check by a direct computation that ker D = Hd (M) . 
We can therefore view II as the spectral projection onto the zero eigenspace of 
D. Consider the unitary group e~ tD , t € R, onl 2 . We shall make use of the 
following essentially well-known formula, 

n= lim - f e- tD dt, (3.3) 

valid in the sense of strong convergence of operators. When checking (|3.3p . let 

2 r 2\ 



u T = ~ [ e - tD dteC(L 2 ,L 
T Jo 



Since 1 1 Hx* 1 1 ,c r_c,2 _t,2 \ < 1, it suffices to check that IIj-x — >■ Ux when x varies in 
a dense subset of I? . We can take this subset to be the set of all finite linear 
combinations of the eigenfunctions of D. To get (|3.3p . we only need to observe 
that when AgR, 

1 f T ux„ fl if A = 0, 



lim - / e UA dt , 
t^+oo T i (^0 ifA^O. 

Now notice that since supp (/) C M_, we have 

e - tD oj f (0) = oj f (t), t > 0, 



and therefore, 



i r T 

IIu/(0)= lim - / <J{t)db. 



We get 



(U k ^ k (0),^ k (0)) L ,= lim 1/ (a/' fe (*),uA fe (0)) 



lim I / (^(0),uA fc (0)) L2 , 



where / t = /(• + t), ft G C^°(M_, C^°(r, AT*M)). Here we have used that 
u/ (t) = u/'(0), as follows from (|2.17|) . An application of Theorem 13 . 1 1 concludes 
the proof. 

□ 



We have the following result with implies Theorem [2? 

Lemma 3.5. Given T C dM and the response operator R^, it is possible to 
construct a finite number of boundary sources fj € Cq°(M_, Cq°(T, AT*M)) 
such that ll k (uh< k (0)) form a basis ofU k D (M), < k < 3. 
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Proof. Let {hj}°° =1 be a dense countable set in Cg°(R_, C^(T, AT*M)). We 
can use the Gram-Schmidt orthogonalization procedure to construct the sources 
fj. More precisely, we define fj recursively by 

to 



3-1 

9i = hj ~ ^Cn fc ^' fc (0), w^CO))^/^ j = 2,3, . . . , 



/: 



(n fe w^' fc (0),u;«' fe (0))y 2 2 ' 



When <7j = 0, we remove the corresponding /ij from the original sequence and 
continue the procedure. The number of sources fj produced by the Gram- 
Schmidt orthogonalization procedure will then be the dimension of 'H^ ) (M), 
according to Corollary 13.31 

□ 



4 Proof of Theorem 12.51 

First notice that as in TheoremEU for any f,h£ C£°(M_, C^(F, A 1 T*M)), the 
knowledge of the response operator Rr allows us to evaluate the inner products, 

(cj f ' k (t),oj h ' k (s)) L 2 , k = 1, 2, t,s>0, 

where u)*,uj are solutions of physical Maxwell's equations f|2. 18|) . 
We have the following controllability result. 

Lemma 4.1. Let uj* be a solution to physical Maxwell's equations (|2.18p . Then 

{n 2 (w /,a (0)) : / € C co (R_,C co (r,A 1 T*M))} = H 2 D (M), 

where H 2 is the orthogonal projection onto the space of the Dirichlet harmonic 
2-fields. 

Proof Let rj 2 G U 2 D (M). Assume that 

(uj f ' 2 (0),r] 2 ) L 2 = for all / e Cg°(R_, C^(T, A 1 T*M)). 
Now (|2.18p and Stokes' formula imply that 

d t (u;f' 2 (t),rj 2 ) L2 = (-^>\tU 2 ) L 2 = -(^(tlS^rfh* - (tu/'^W) 

= -(f(t),n V 2 ). 

Thus, 

(f(t),n V 2 ) = -(^ 2 (0),77 2 ) L2 + (J> 2 {-T f W) L . = 0, 
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for all / <G C^°(M_, C °°(r, A^M)). Hence, nr/ 2 = on T. Moreover, Arj 2 = 
on M and trj 2 = on dM. By the unique continuation principle, we get rf = 0. 

□ 

Proceeding further as in Subsection 13.3^ we can recover the first absolute Betti 
number (3\(M) from the knowledge of V and Ry- This completes the proof of 
Theorem 12.51 
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